Vortex States Induced by Proximity Effect in Hybrid 
Ferromagnet-Superconductor Systems 
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We consider superconductivity nucleation in multiply connected mesoscopic samples such as thin- 
walled cylinders or rings placed in electrical contact with a ferromagnet. The superconducting 
critical temperature and order parameter structure are studied on the basis of linearized Usadel 
equations. We suggest a mechanism of switching between the superconducting states with different 
vorticities caused by the exchange field and associated with the oscillatory behavior of the Cooper 
pair wave function in a ferromagnet. 
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I. INTRODUCTION 



The origin of the vortex states in hybrid ferromag- 
net (F)-supcrconductor (S) structures is closely related 
to the basic mechanisms responsible for the interplay 
between the ferromagnetic and superconducting order- 
ings (see, e.g., Refi). The first mechanism is associated 
with the orbital effect, i.e., interaction of Cooper pairs 
with the magnetic field induced by magnetic moments 2 -. 
This field causes the appearance of inhomogeneous super- 
conducting phase distributions and spontaneous vortex 
states. The switching between these states characterized 
by different winding numbers can result in an oscillatory 
behavior of the critical temperature T c as a function of 
the external field H, which resembles the Little-Parks 
effect in multiply connected superconducting samples^. 
Such nonmonotonic behavior of T C {H) was shown to be 
inherent to hybrid F/S systems with magnetic dots or 
domains^^, that create a "magnetic template" for nu- 
cleation of superconducting order parameter. The second 
mechanism arises from the exchange interaction which 
comes into play because of the proximity effect, when 
the Cooper pairs penetrate into the F layer and induce 
superconductivity there. The latter effect is known to 
result in the damped-oscillatory behavior of the Cooper 
pair wave function in a ferromagnet^ which is the cause 
of a number of fascinating interference phenomena in hy- 
brid F/S structures. In particular, this peculiar proxim- 
ity effect reveals itself in the oscillating^ or re-entrant^ 
behavior of the critical temperature as a function of a fer- 
romagnetic layer thickness in layered F/S structures, and 
is responsible for the formation of Tr-junction o 10 ^ 1 ' 12 . A 
Joscphson 7r-junction is a generic example of the system 
where the proximity effect in a ferromagnetic subsystem 
is used to obtain an energetically favorable superconduct- 
ing state with a nontrivial distribution of the order pa- 
rameter phase. A resulting distinctive feature of the sys- 
tems with 7r-junctions is a possible unusual ground state 
with spontaneous supercurrents and, in particular, with 
spontaneously formed vortices. For a Joscphson junc- 
tion with a step-like change in the F layer thickness such 



spontaneously formed vortex states have been discussed, 
e.g., in RefsJ&Iiils. 

It is the purpose of this paper to examine a possibility 
to realize a switching between the spontaneously created 
vortex states in multiply connected samples caused by 
the proximity effect with a ferromagnet. We focus on 
the behavior of critical temperatures for superconducting 
states with different vorticities and, thus, in some sense 
study an analog of the Little-Parks effect caused by the 
exchange interaction mechanism. 

To elucidate our main results we start from a qual- 
itative discussion of the proximity effect on the super- 
conducting ground state in a thin-walled superconduct- 
ing shell surrounding a cylinder (core) of a ferromagnetic 
metal (see Fig. [TJ. We expect that the superconduct- 
ing ground state in such geometry should be strongly 
influenced by the damped-oscillatory behavior of the su- 
perconducting order parameter in ferromagnetic cylinder 
and, thus, controlled by the ratio of the period of the or- 
der parameter oscillations (~ £/) to the radius Rf of the 
F core. Indeed, for Rf < £/, the variation of the pair 




FIG. 1: The schematic representation of the F/S system un- 
der consideration: thin-walled superconducting shell around 
a ferromagnetic cylinder. Here Rf is the radius of the F core, 
and R a is the outer radius of the S shell, (r, 0, z) is the cylin- 
drical coordinate system. 
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FIG. 2: The schematic behavior of the superconducting order 
parameter inside the F cylinder, (a) The curve ^{x) repre- 
sents sketchy the behavior of the pair wave function in the 
L — phase. Due to symmetry the derivative d x ^ is zero at 
the center of F cylinder, (b) The pair wave function in the 
phase with L = 1 vanishes at the center of F cylinder, and 
ty(x) has a 7r— shift in diametrically opposite points. 



II. MODEL 

We assume the elastic electron-scattering time r to be 
rather small, so that the critical temperature T c and ex- 
change field h satisfy the dirty-limit conditions T c t <C 1 
and ftr < 1. In this case a most natural approach to 
calculate T c is based on the Usadel equations^ for the 
averaged anomalous Green's function Ff and F s for the 
F and S regions, respectively (see^ for details). Near the 
second-order superconducting phase transition, the Us- 
adel equations can be linearized with respect to the pair 
potential A(r). In the F (S) region these linearized Us- 
adel equations take the form 

-^V 2 F f + (\u\+ih S gnu})F f =0, (1) 



-yV 2 F s + \lj\F s = A(r) 



(2) 



wave function *F along a line crossing the F core appears 
to be modest and the order parameter can not change its 
sign along the line. It means that in the ground state the 
superconducting phase in diametrically opposite points 
must be the same. The resulting angular momentum L 
of the pair wave function 'F is equal to zero (see Fig. [2K)- 
This state is analogous to the 0— phase state of SFS lay- 
ered structures^. For a larger radius Rf > £/, the pair 
wave function may cross zero at the axis of the F cylin- 
der which causes a it— shift in the superconducting phase 
in diametrically opposite points. In this case the angu- 
lar momentum L of the pair wave function ^ is nonzero 
(Fig. [2)3). Thus, the penetration of Cooper pairs into the 
F core and exchange interaction can induce the super- 
conducting states with a nonzero vorticity. It is natural 
to expect that the supercurrents flowing in such states 
with L / are responsible for an additional contribu- 
tion to the vortex energy. The interplay between the 
exchange effect and supercurrent depending energy term 
may result in a subsequent switching between the states 
with different vorticities, as the F core radius increases. 
An obvious consequence of these transitions between the 
states with different L should be a nonmonotonic depen- 
dence of the critical temperature T c on the F core radius 
and exchange field. 

The paper is organized as follows. In Sec. [IT] we briefly 
discuss the basic equations. In Sec. |TTI] we study the 
switching between different vortex states in two model 
F/S systems. The first system consists of a thin- walled 
superconducting cylindrical shell surrounding a cylinder 
of a ferromagnetic metal. The second one is a planar 
structure which consists of a mcsoscopic superconducting 
ring placed at the surface of a thin ferromagnetic film. 
For both cases we assume that there is a good electrical 
contact between the F and S regions, to assure a rather 
strong proximity effect. We summarize our results in 
Sec. ED 



The superconducting critical temperature T c is deter- 
mined from the self-consistency condition for the gap 
function: 

A(r)^ + ,T c ^(^--F sM )=0. (3) 

Here Df and D s are the diffusion constants in the fer- 
romagnet and superconductor, respectively, and uj = 
(2n + 1)ttT c is a Matsubara frequency at the tempera- 
ture T c . Equations (UJ),© must be supplemented with 
the boundary condition at the outer surfaces 

d*F Ls = , (4) 

and at the interface between the F and S metals>i£ 

a s d n F s = <j f d n F f ; F s = F f + 7fc 6 d n F f . (5) 



Here £ s = ^D s /2ttT c0 is the superconducting coherence 
length, tjf and a s are the normal-state conductivities of 
the F and S metals, 71, is related to the F/S boundary 
resistance Rt per unit area through 7b£ s = Rb<Jf, and d n 
denotes a derivative taken in the direction perpendicular 
to the outer surfaces or to the F/S interface. For the 
sake of simplicity we assume h 3> irT c o and neglect the 
proximity effect suppression caused by a finite F/S inter- 
face resistance^, i.e., take 7;, — > 0. In this regime we get 
Ff = F s at the F/S interface. 

For a system with a cylindrical symmetry the vorticity 
parameter L just coincides with the angular momentum 
of the Cooper pair wave function. Choosing cylindrical 
coordinates (r, 0, z) we look for solutions of the equations 
(H]),©,© characterized by certain angular momenta L: 

A(r)=A(r,z)e* Le , f>,.(r) = /,>, z) e M , (6) 

According to the equations ©,©,© there is a symme- 
try Ff >s (ui) = Ft 8 (—u>), so that we can treat only positive 
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lo values. The Usadcl equations (HJ) , ([2j) can be written in 
the form 

+ihf f = 0, 
D x (\ „ , „ _ „„ „ L 2 

,•2 



2 \r 



9rMr/ S ) + 9 z 2 / S - 3 / S ) 

+uj f s = A . 



An appropriate self-consistency equation (|3|) can be 
rewritten as follows: 

Aln^+27rT c V (--Ref s (u)) = 0. (9) 




III. CRITICAL TEMPERATURE OF VORTEX 
STATES 

Now we proceed with the critical temperature calcu- 
lations for different vortex states. For the sake of defi- 
niteness we consider here two generic examples of hybrid 
F/S systems which we believe to manifest the vorticity 
switching scenario suggested in the introduction. 



A. Thin-walled superconducting shell around a 
ferromagnetic cylinder 

Consider a superconducting cylindrical shell of a thick- 
ness W = R s — Rf <C Rf surrounding a thin cylinder 
of a ferromagnetic metal with a uniform magnetization 
M = Mz n . Here Rf is the radius of the F core, and R s is 
the outer radius of the S shell (see Fig. [1}. Naturally, to 
observe the pronounced influence of the proximity effect 
on the transition temperature, the thickness of the S shell 
W must be smaller than the superconducting coherence 
length £ s . 

For a thin and long F cylinder with the magnetization 
direction chosen along the z— axis the magnetic field B 
vanishes outside the F region. However, this magneti- 
zation induces a vector potential Ag ~ 2nMRf in the S 
region, which may result in a standard Little-Parks ef- 
fect of electromagnetic origin. We may take account of 
this vector potential in Eq. ([5]) replacing the vorticity pa- 
rameter L by the value L — <3>/$o, where 3>o = 27rSc/e 
is the magnetic flux quantum and <3? = 4ir 2 MRj is the 
total magnetic flux. For the standard Little-Parks effect 
the critical temperature vs $ oscillates with a period $o 
and an amplitude AT C ~ T c £ 2 /i? 2 . For typical param- 
eters M ~ 10 2 G, T ~ 10 K, D s ~ 10cm 2 /s and R f of 
order of several £/ ~ 10 nm lengths we get $/$o "C 1 
and AT C <C T c . These simple estimates allow us to find 
a region of parameters where we can exclude the effect 
of magnetic field on T c . The above effect of the mag- 
netic field can be also weakened provided we decrease the 
height of the F/S cylinder going over to the case of a thin 



FIG. 3: The dependence of the critical temperature T c on the 
F core radius Rf for two values of the vorticity L — (solid 
line) and L = 1 (dashed line). Here we choose W = 0.5£ s ; 
a s /a f = 2.5; = 0.265. 



disk when the field B is suppressed due to the demagne- 
tization factor. Note that choosing the magnetization di- 
rection in the plane perpendicular to the cylinder axis we 
can get rid of the standard Little-Parks effect completely 
because of the absence of the magnetic field component 
along the cylinder axis. On the contrary, the vorticity 
switching scenario studied below does not depend on the 
magnetization direction. 

We look for homogeneous along z solutions of the equa- 
tions ([7]), ijHJ) with a certain angular momentum L. In 
this case the equation ([7|) in the F cylinder can be read- 
ily solved: 



ff = CI L (q f r), q f 







(10) 



Here Il{u) is the modified Bessel function of first kind 
of order L, and £/ = y/Df/h is the characteristic length 
scale of the order parameter variation in the F metal. In 
the dirty limit parameter £f determines both the length 
scale of oscillations and decay length for the Cooper pair 
wave function in a ferromagneti. The boundary condi- 
tions ©,© for Eq. © take the form: 



s dr 
a L = 



L 



a L q f (Tff s (Rf) , ^ 



0, (11) 



u 



f 



> u f = if R f- 



For W <C £ s , the variations of the functions / s (r) and 
A(r) in the superconducting shell are small: f s (r) — f, 
A(r) ~ A. Therefore, we can average Eq. © over the 
thickness of the S shell, using the boundary conditions 
(fTTj) to integrate the term d r (rd r f s ). Finally, we obtain 
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FIG. 4: The typical dependences of the critical temperature T c on the F core radius Rf for different values of the vorticity 
L = (solid line), L = 1 (dashed line) and , L — 2 (dotted line). Here we choose W = 0.5£ s ; = 0.28 and different values 

of the ratio a a jaf. a) a s /of = 3; b) cr s /<x/ = 2.7; c) a s /af = 2.5; d) o B joj = 2. 



the following expression: 



/ = 



w ■ 



^7 



<* ill 
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(12) 



Substituting Eq. (fT2")) into Eq.© one obtains a self- 
consistency equation for the critical temperature T c : 



In — = # | - ) - Re # f - + fl L 

Too \2j \2 



(13) 



where \& is the digamma function. The depairing param- 
eter 



_ 1 T c0 2 



^7 



OIL 



(14) 



is responsible for the superconductivity destruction in the 
shell due to both the exchange effect and the supercurrent 
flowing around the cylinder. 



Figure [3] shows a typical dependency of the critical 
temperature T c on the F core radius Rf, obtained from 
Eqs. (fT3|) . (fT4|) for different winding numbers L. We see 
that for a small F cylinder radius fl/ <C only the state 
with L — appears to be energetically favorable. The 
influence of the proximity effect is weak and the criti- 
cal temperature T c is close to T c $. The T c of a vortex 
state with L ^ is suppressed because of a large su- 
percurrent energy. The increase in the radius Rf results 
in a decrease in T c for the state with L = and re- 
duce the kinetic energy of supcrcurrents for L ^ 0. At 
the same time, the damped-oscillatory behavior of the 
superconducting order parameter in a fcrromagnet be- 
comes important. If the diameter of the F cylinder is 
comparable with the period of the order parameter oscil- 
lations (~ £/) and satisfies approximately the conditions 
of the 7T— phase superconductivity in layered F/S struc- 
tures (2£f < 2Rf < 5£/)i, then there appears a 7r— shift 
in the phase of the superconducting order parameter in 
diametrically opposite points (see Fig. (2b) . In this case, 
the critical temperature of L = 1 state becomes higher 
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FIG. 5: Superconducting ring lying at the surface of a thin 
ferromagnetic film. Here Ro and W are the radius and the 
width of the S ring, and df (d a ) is the thickness of a ferro- 
magnetic (superconducting) layer. 



than the critical temperature of the state with L = 0. 
Thus, our calculations confirm the qualitative arguments 
given in the introduction. The penetration of Cooper 
pairs into the F core and the phase shift of the pair wave 
function \I/ due to the exchange interaction can induce 
vortex states in the superconducting shell. 

To illustrate the scenario of switching between the 
states with different vorticities L and a nonmonotonic de- 
pendence of the critical temperature T c vsRf we present 
here several T c (Rf) curves for various conductivity ratio 
(see Fig. [4j. We see that the ratio o s /<Jf of the normal- 
state conductivities of the F and S metals is an important 
factor, controlling the generation of the vortex states in 
the F/S structure under consideration. 



B. Superconducting ring at the surface of a thin 
ferromagnetic film 

As a second example we consider a planar hybrid sys- 
tem, i.e., a superconducting ring lying at a ferromag- 
netic film with a uniform in-plane magnetization M (sec 
Fig. [S]). Such version of the setup can be more con- 
venient for the experimental observation of the switch- 
ing phenomena. The superconducting ring of the ra- 
dius Ro and the width W <C Ro occupies the region 
Ro - W/2 < r < Ro + W/2, < z < d s . The 
boundary conditions at the interfaces with vacuum yield 
d r f s (R ± W/2, z) = 0, and d z f s (r,d 8 ) = 0. For sim- 
plicity, we consider only the case of a rather thin ring 
with W <C £/ lS and d s <C £ s which allows us to as- 
sume the variations of the functions / s (t") and A(r) in 
the superconducting ring to be small. Thus, we can av- 
erage Eq. ([8]) over the volume of the S ring, integrate the 
terms d r {rd r f s ) and 9f/ s , and make use of the bound- 
ary condition at the interface with vacuum. Finally, we 
get the following expression for the derivative d z f s at 
r = Ro, z = 0: 



±d z f.= A(A- W /.)- 

d s D s V Ro 



u 



(15) 



The F metal occupies the region —df < z < 0. We 
will address the case of a very thin F film: df <C £/. The 
boundary conditions at the interfaces with vacuum yield 
d z ff = 0. Far from the superconducting ring, i.e. for 
r - Ro - W/2 > df and for R - W/2 - r > df we can 
average the Usadel equation ([7|) over the thickness df. 



D I 



d r (rd r ff) 



+ ihf f = Q. (16) 



The solution of this equation reads: 

CxI L (q f r) , Ro-W/2-r^df 



ff(r) 



C 2 K L {q f r) 



Ro - W/2 > df 



(17) 



where Il{u) and K~l(u) are the modified Bessel functions 
of order L. Making use of this solution one can easily get 
the following relations between the function ff and its 
derivative d r ff at r = Rq ± e: 



dr 
dr 



R -e 



Ro+e 



If 

Qf 



L 
u 
L 
u 



I 



L + l 



(uq) 



Il{uo) 
K l+1 (u q ) 
K L {uo) 



ff(Ro-e), (18) 
f f (Ro+s), (19) 



where u = q/Ro and max[Fy, df] <£<^j< i?o- As- 
suming ff to be a slow function of r one can write a 
boundary condition on the derivative jump: 



dr 



Ro+e 



Ro-e 



Q L = 



-qfQLff(Ro), 

Il+i(u ) . K L+1 (uo) 



(20) 



Il{u ) 



K L {uo) 



On the other hand in the region \r — Rq\ <C £/ the Eq. ([7]) 
takes a simple form: 

1 -d r {rd r ff) + d 2 Jf = Q. 

Integrating this equation over the region Rq — e < r < 
Ro + s and over the ferromagnetic film thickness and 
making use of the boundary conditions described above 
we obtain: 



Wd z f f 



Qfdf QLff(Ro) 



(21) 



z=0,r=R 



As before we restrict ourselves to the case of low F/S in- 
terface resistance assuming 76 = in ([5|). In this regime, 
f f {R o ,0) ^ fs[Ro,0) = f. The Eqs. {T®,^ and the 
boundary conditions at the F/S interface ([5]) determine 
the amplitude /: 



/ 



CO 



Ds 
2 



L 
Ro 



If 
rjW 



(22) 
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FIG. 6: The dependence of the critical temperature T c on the S ring radius Ro for different values of the vorticity L = (solid 
line), L = 1 (dashed line) and , L = 2 (dotted line). Here we choose d s /df = 1, W = 0.5£ s , = 0.1: a) <J 3 /<jf = 2.5; b) 

a s /a f = 2.1; c) a B /a f = 2.03; d) a B /a f = 2.0. 



where = a s d s /(Tfdf. Substitution of (f2"2"| into Eq.@ 
results in the self-consistency equation for the critical 
temperature T c of the F/S hybrid (fT3|) . where the de- 
pairing parameter is determined by the following ex- 
pression: 



o 1 T °° c 2 



L 



If 
Wr) 



(23) 



In Fig.[6j we present typical dependences of the critical 
temperature T c on the S ring radius Rq for the different 
orbital number L, obtained from Eqs. (fl~3|) . (|23|) . The 
curves appear to be qualitatively similar to the ones ob- 
tained for a superconducting thin-walled cylinder in the 
previous subsection. Note that for a particular choice of 
parameters (see Fig|6j:) T c vanishes in a certain interval 
of i?o values and we observe an interesting re-entrant be- 
havior of the critical temperature for a state with zero 
vorticity. 



IV. SUMMARY 



To summarize, we suggest a mechanism of switching 
between the superconducting states with different vor- 
ticities caused by the exchange field in the hybrid S/F 
structures and associated with the oscillatory behavior 
of the Cooper pair wave function in a ferromagnet. We 
defined the range of the system parameters at which the 
predicted effect can be experimentally observable. The 
most restrictive condition is imposed on the relation be- 
tween the superconducting cylinder or ring radius Rq and 
coherence lengths £ s and £/. On the one hand the radius 
should essentially exceed the superconducting coherence 
length to decrease the kinetic energy of supercurrents in 
the states with nonzero vorticity, but on the other hand 
it should be of the order of only several £/ lengths to 
ensure a rather strong influence of the proximity effect. 
As a result, to observe the switching effect we need to 
consider the S/F systems with a rather large ratio £/ /£ s . 
This ratio can be increased if we choose a superconduct- 
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ing material with a rather short coherence length, e.g., 
heavy fermion compounds^. Despite of such restriction 
we believe that the proximity induced switching between 
the vortex states can be experimentally observable and 
would provide an interesting manifestation of interference 
effects in mesoscopic superconductivity. 
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